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We present the results for a chemical abundance analysis between planet-hosting and stars without planets for 12 refractory elements 
for a total of 1 1 1 1 nearby FGK dwarf stars observed within the context of the HARPS GTO programs. Of these stars, 109 are known to 
harbour high-mass planetary companions and 26 stars are hosting exclusively Neptunians and super-Earths. We found that the [X/Fe] 
ratios for Mg, Al, Si, Sc, and Ti both for giant and low-mass planet hosts are systematically higher than those of comparison stars at 
low metallicities ([Fe/H] < from -0.2 to 0. 1 dex depending on the element). The most evident discrepancy between planet-hosting and 
stars without planets is observed for Mg. Our data suggest that the planet incidence is greater among the thick disk population than 
among the thin disk for mettallicities bellow -0.3 dex. After examining the [a/Fe] trends of the planet host and non-host samples we 
conclude that a certain chemical composition, and not the Galactic birth place of the stars, is the determinating factor for that. The 
inspection of the Galactic orbital parameters and kinematics of the planet-hosting stars shows that Neptunian hosts tend to belong to 
the "thicker" disk compared to their high-mass planet-hosting counterparts. We also found that Neptunian hosts follow the distribution 
of high-or stars in the UW vs V velocities space, but they are more enhanced in Mg than high-o- stars without planetary companions. 
Our results indicate that some metals other than iron may also have an important contribution to planet formation if the amount of 
iron is low. These results may provide strong constraints for the models of planet formation, especially for planets with low mass. 

Key words, stars: abundances - stars: planetary systems 



1. Introduction 

Since the first discovery of an exoplanet around a solar-like star, 
(Mayor & Queloz 1995), more than 750 exoplanets have been 
discovered. More than one hundred of these planets have masses 
less than 30 M @ . With the ever increasing number comes the real- 
ization that the nature of these planets is diverse, thereby adding 
complexities to an already complicated problem. 

There are essentially two main and competitive planet forma- 
tion models: the core-accretion model (e.g. Pollack et al. 1996; 
Ida & Lin 2004; Mordasini et al. 2009) and the disk instability 
model (e.g. Boss 1997;2002). The tendency for planets to orbit 
metal-rich stars (Gonzalez 1998; Gonzalez et al. 2001; Santos 
et al. 2001, 2003, 2004; Fischer & Valenti 2005; Sousa et al. 
2008; Neves et al. 2009; Johnson et al. 2010; Sousa et al. 201 la; 
Adibekyan et al. 2012 submitted) lends strong support to the 
first scenario, according to which planets grow through accretion 
of solid, metal-rich material to form massive cores. The discov- 
ery of several planets orbiting metal-poor stars (e.g. Cochran et 
al. 2007; Santos et al. 2010) shows, however, that giant planet 
formation is not completely inhibited in the metal-poor regime 
(see also discussion in Santos et al. 2004, 2011). Recent studies 
showed that the well-established metallicity correlation seems to 
be only observed for giant planets. The stars hosting Neptunes 



* Based on observations collected at the La Silla Parana Observatory, 
ESO (Chile) with the HARPS spectrograph at the 3.6-m telescope (ESO 
runs ID 72.C-0488, 082.C-0212, and 085.C-0063). 



and super-Earth-class planets show typical metallicities of the 
stars without planets (e.g. Udry et al. 2006; Sousa et al. 2008; 
Ghezzi et al. 2010; Mayor et al. 2011; Sousa et al. 2011a). This 
gives us interesting hints about the planet formation processes. 

In the conventional core-accretion scenario, the formation of 
planetesimals starts from the condensation of heavy elements 
(metals). In this context, the study of the heavy element abun- 
dances is very important. Previous studies that aimed to clarify 
whether the planet-hosting stars (PHS) are different from stars 
without planets in their content of individual heavy elements 
(other than iron) yielded clontradictory results. Although most 
studies showed no significant differences in the overall trends of 
[X/Fe] between PHSs and stars without any known planetary- 
mass companions (e.g. Takeda 2007; Bond et al. 2008; Neves et 
al. 2009; Gonzalez Hernandez et al. 2010; Delgado Mena et al. 
2010), some works have reported possible enrichment in some 
species (Gonzalez et al. 2001; Santos et al. 2000; Sadakane et 
al. 2002; Bodaghee et al. 2003; Fischer& Valenti 2005; Beirao 
et al. 2005; Gilli et al. 2006; Bond et al. 2006; Robinson et al. 
2006; Gonzalez & Laws 2007; Brugamyer et al. 2011; Kang et 
al. 2011). 

Although theoretical modeling suggests that metallicity is 
a key parameter of planet formation, Haywood (2008), study- 
ing the memberships of PHSs to different stellar populations, 
proposed that the presence of giant planets might be primarily 
a function of a parameter linked to galactocentric radius, but 
not metallicity and the apparent correlation between metallic - 
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ity and the detection of planets is a natural consequence of that. 
Haywood (2009) suggested that this parameter, which depends 
on the distance from the Galactic Center, could be the density of 
molecular hydrogen. If the rate of giant planets does not depend 
on metallicity, then the core-accretion theory of planet formation 
(e.g. Ida & Lin 2005; Mordasini et al. 2012) loses its most im- 
portant observational support, and in turn strengthens the gravi- 
tational instability theory (e.g. Boss 2001; Cai et al. 2006). The 
kinematics of extra-solar planet hosts and their relation to differ- 
ent stellar populations and moving groups have been discussed 
in some other works (e.g. Gonzalez 1999; Reid 2002; Barbieri & 
Gratton 2002; Santos et al. 2003; Ecuvillon et al. 2007; Neves et 
al. 2009; Gonzalez 2009). The results do not allow one to reach 
any clear conclusion. 

The study of extrasolar planets requires very high quality 
data. In particular, very high-precision radial-velocity measure- 
ments are needed to detect planets. Likewise, finding possible 
abundance differences (sometimes a very subtle task) between 
stars with and without planets also requires using large stellar 
samples with accurate and homogeneous abundance determina- 
tions. In this paper, we present a uniform spectroscopic analysis 
of 1 1 1 1 FGK dwarfs observed within the context of the HARPS 
GTO planet search program. Of these stars, 109 are known to 
harbour high-mass planetary companions and 26 stars are host- 
ing exclusively Neptunians and super-Earths. The large size and 
homogeneity of this sample enables a very robust differential 
comparison of the chemical abundances of the stars with and 
without planets with minimal internal uncertainties. The rela- 
tively high number of low-mass planet hosts in the sample al- 
lowed us to study their chemical properties separately from their 
Jupiter-mass planet hosting counterparts. The introduction of the 
sample and the methods of the chemical abundance determina- 
tion and analysis are described in details in our previous paper 
(Adibekyan et al. 2012, submitted). This paper is organized as 
follows: In Sect. 2 we briefly introduce the physical properties 
of the stars in the sample. The study of the abundances of re- 
fractory elements relative to iron in exoplanet-hosting stars can 
be found in Sect. 3. The kinematical properties of the exoplanet- 
hosting stars is presented in Sect. 4. Finally, in Sect. 5 , we draw 
our main conclusions. 



MOOG 1 (Sneden 1973) and a grid of Kurucz ATLAS9 plane- 
parallel model atmospheres (Kurucz et al. 1993). The reference 
abundances used in the abundance analysis were taken from 
Anders & Grevesse (1989). The line list and atomic parameters 
of Neves et al. (2009) were used adding one Cai line and ex- 
cluding five Nil lines, two Sil lines, two Till lines, and five Til 
lines (Adibekyan et al. 2012, submitted). The equivalent widths 
were automatically measured with the ARES 2 code (Automatic 
Routine for line Equivalent widths in stellar Spectra - Sousa et 
al. 2007). The final abundance for each star and element was 
calculated to be the average value of the abundances given by all 
lines detected in a given star and element. Individual lines for a 
given star and element with a line dispersion higher than a fac- 
tor of two than the rms were excluded. The total uncertainties 
in the [X/H] abundances and [X/Fe] ratios do not exceed 0.04 
and 0.03 dex for stars with r e jf = r o +500 K, respectively, and 
are less than 0. 1 dex (except Scl, Til and VI) for the stars with 
temperatures very different from that of the Sun. 

With our large stellar sample, we were able to detect and cor- 
rect the [X/Fe] trends with T e e for some elements. For All and 
Col we observed a systematic trend with T e g in all temperature 
ranges, and for Til, Scl, V, CrII, and Na we observed a trend with 
r e ff in the low-temperature domain. 

As a check of our method and analysis, we compared our 
derived abundances with those obtained by Bensby et al. (2005), 
Valenti & Fisher (2005), Gilli et al. (2006), and Takeda (2007) 
for stars in common with these works. In genereal we found a 
good agreement with these previous studies. 

For more details about the sample, abundance determination 
methods, analysis, and errors we refer the reader to Adibekyan 
et al. (2012, submitted). 



3. Abundances in planet hosts 

In this section we perform a detailed analysis of the [X/Fe] 
distributions of stars hosting Jupiter-mass planets, stars with 
Neptunes and super-Earths, and stars without planets. At the 
end of this section we investigate the connection of metal-poor 
planet-hosting stars to the Galactic thick disc. 



2. The sample 

The sample used in this work consists of 1111 FGK stars ob- 
served with the HARPS spectrograph (Mayor et al. 2003) at the 
ESO 3.6 m telescope (La Silla, Chile). The stars are slowly ro- 
tating, non-evolved, and in general have a low level of activity. 
Most stars have effective temperatures 4600 K < T e g < 6300 K 
(only 4% are outside of this region) surface gravities 4 < log g < 
5 dex (the number of "outliers" again is very small, only 25 stars 
with log g < 4 dex) and they lie in the metallicity range of -1.39 
< [Fe/H] < 0.55 dex. 

Precise stellar parameters for all stars were determined in the 
same manner and from the same spectra as used in our previous 
study. For details we refer the reader to Sousa et al. (2008, 201 la, 
201 lb). The typical uncertainties in the atmospheric parameters 
are of about 30 K for T & s , 0.06 dex for log g, and 0.03 dex for 
[Fe/H]. 

Elemental abundances for 12 elements (Na, Mg, Al, Si, Ca, 
Ti, Cr, Ni, Co, Sc, Mn, and V) were determined using a lo- 
cal thermodynamic equilibrium (LTE) analysis relative to the 
Sun with the 2010 revised version of the spectral synthesis code 



3. 1 . [X/Fe] of planet host stars 

Through studing the [X/Fe] histograms of exoplanet-hosting 
starsone can characterize the distribution of the individual el- 
ements relative to iron for the sample in general, and find the 
limits of the distributions. The [X/Fe] distributions and the cor- 
responding cumulative fractions for the total sample are depicted 
in Fig. 1 and Fig. 2, respectively. In the plots (and in the sub- 
sequent figures) we used the average of Til & Till for Ti, the 
average of CrI & CrII for Cr, and the average of Scl & Sell for 
Sc to increase the statistics. The abundance ratios do not change 
when using the mean values as compared to using the different 
ions separately. 

The inspection of Fig. 1 and Fig. 2 suggests that, in gen- 
eral, the [X/Fe] distributions of Neptunian hosts (NH) start at 
higher [X/Fe] values compared to the distributions of their gi- 
ant planet-hosting counterparts. For most elements this shift in 
part reflects the fact that Neptunian and Jovian hosts (JH) have 



1 The source code of MOOG2010 can be downloaded at 
http://www.as.utexas.edu/~chris/moog.html 

2 The ARES code can be downloaded at 
http://www.astro.up.pt/sousasag/ares 
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Fig. 1. [X/Fe] distribution of the different elements. The stars with giant planets and without planets are represented by a red dashed and black dot- 
ted lines, respectively. The stars that exclusively host Neptunians and super-Earth planets are represented by a shaded blue. The Neptunian/super- 
Earth and Jovian hosts distributions were set smaller (0.5 and 0.8 times, respectively) for the sake of clarity. 



a different metallicity distribution and hence different [X/Fe] ra- 
tios because of the Galactic chemical evolution trends. However, 
this is not the case of Mg, for which the shift of the distribution 
toward the higher [X/Fe] values is most evident. This shift exists 
even when we compare the two subsamples in the same metallic- 
ity interval. The Kolmogorov-Smirnov (K-S) statistics predict a 
10~ 4 probability (P) that Neptunian and Jovian host stars have 
the same [Mg/Fe] distribution. The K-S tests also give signifi- 
cant probabilities that the [X/Fe] distribution of the two planet- 
hosting samples are different for Ca (K-S P ~ 10~ 4 ) , Ti (K-S P 
* 0.04), and Mn (K-S P ~ 0.05). 

Table 1 lists the average values of [X/Fe] for three groups of 
stars, along with their rms dispersion, the number of stars used 
in their determination, and the difference of averages between 
Neptunian and Jovian hosts and stars without planets. These dif- 
ferences range from about -0.005 (Cr) to 0.052 (Mg) for NHs. 
The maximum overabundance of JHs are found for manganese 
A[Mn/Fe] « 0.06 dex. We note that previous studies (Bodaghee 
et al. 2003; Zhao et al. 2002) have also suggested the differ- 
ence in [Mn/Fe] between PHSs and comparison stars, though 
the difference was statistically within their scatter. As mentioned 
above, Kang et al. (201 1) also claimed to have found even signif- 
icant overabundances of [Mn/Fe] of planet-hosting stars with the 



K-S probability 0.0015% that the host and non-host stars have 
the same [Mn/Fe] distribution. We note that these differences can 
be partially "affected" by the underlying [Fe/H] distributions of 
planets with different mass and control sample. 

To evaluate the probability that the abundances of the two 
planet-hosting groups and the sample without planets have the 
same [X/Fe] distribution, we performed a K-S test. To avoid 
confusion from the multiple values of [X/Fe] (for "a-like" el- 
ements) at low metallicities due to the thin and thick disk (e.g. 
Bensby et al 2003; Adibekyan et al. 201 1) and to remove (at least 
partially) the differences of the underlying iron distributions for 
the three samples, we established a cutoff in [Fe/H]. Only stars 
with metallicities between -0.2 and +0.4 dex were used for the 
K-S tests. The values of the probability that the [X/Fe] distri- 
butions belong to the same population (Prob) are presented in 
the Table 2. We note that, in general, the Prob values did not 
change when we extended the metallicity region from +0.4 to 
+0.55 dex (+0.55 dex is the metallicity limit for our sample), 
but they decreased (especially for "a-like" elements) if we added 
the stars with metallicities from -0.2 to -0.6 dex (-0.6 dex is the 
lower metallicity limit for planet hosts in our sample). This is 
expected because in the low-metallicity regime almost all planet 
hosts have high [a/Fe] values, typical for thick-disk stars, while 
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Table 1. Average abundance ratios [X/Fe] for stars without planets, with giant planets, and stars that exclusively host Neptunians, along with 
their rms dispersion, the number of stars used in their determination, and the difference of averages between Neptunian and Jovian hosts and stars 
without planets. 



Species 


Jovian hosts 




Neptunian hosts 




Without planets 




Difference of averages 
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<[X/Fe]> 


<r 


N 


<[X/Fe]> 


cr 


yv 


<[X/Fe]> 


cr 


N 


Jovian - 
Non-hosts 


Neptunian - 
Non-hosts 


Na 


0.054 


0.076 


109 


0.061 


0.069 


26 


0.031 


0.098 


975 


0.023 


0.030 


Mg 
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0.083 
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0.016 
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Cr 
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976 


-0.017 
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Mn 
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Fig. 2. Cumulative distribution functions of [X/Fe]. The stars with 
giant- and low-mass planets are represented by a red dashed and a blue 
solid line, respectively. The stars without planets are represented by a 
black dotted line. 



stars without planets belong to the thin and thick disks with low 
and high [a/Fe] values, respectively (see Sect 3.2 and 3.3). For 
the metallicity distribution, the probability that the [Fe/H] distri- 
butions of NHs and non-hosts belong to the same population is 
about 99% and for JHs and non-hosts the test gives 0% proba- 
bility. The K-S probabilities for NHs reflects the real differences 
in [X/Fe] between these and non-host stars because their under- 
lying iron distributions fully match. The K-S test implies that 
the stars hosting Neptunians and stars without planets are in- 
deed from separate groups in terms of [Mg/Fe], [Si/Fe], [Sc/Fe] 
and [Al/Fe]. For Jupiter-mass planet hosts the K-S test implies 
that they are different from stars without planets by their [X/Fe] 
distribution for all the elements (except Ti and V) although we 
should note that their [Fe/H] distributions are different as well. 

Fig. 3 illustrates the [X/Fe] distribution from a different per- 
spective: the histograms of the number of stars with giant (red 
dashed) and low-mass planets (shaded blue) compared to the to- 
tal number of stars of each bin (0.05 dex). We note that some 
bins with extreme [X/Fe] values show an unreal 100% percent- 
age with only one star in the bin. For the sake of clarity we cut 
the y axis at 30%. For Co, Na, Ni, V, and Mn we observe that 



o 
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[X/Fe] 

Fig. 3. Percentage of stars with giant (red dashed) and exclusively 
Neptunian and super-Earth (shaded blue) planets as a function of 
[X/Fe]. For the sake of clarity we cut the y axis at 30 %. 



there is a general increase in the percentage of stars with gi- 
ant planets, with increasing [X/Fe]. An opposite trend shows the 
distributions of Ca and Cr. It is difficult to claim that there is 
an opposite trend for Mg because the distribution is in general 
flat when we exclude the lowest [Mg/Fe] bin with the "unreal" 
100% frequency. We also found that the probability to find a 
Neptunian mass planet is increasing if the host has high [X/Fe] 
values for all studied elements except Mn. Interestingly, most of 
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Fig. 4. [X/Fe] abundance ratios against [Fe/H] for the stars with and without planets. The symbols and error bars indicate the average and standard 
deviation, respectively, of each bin (0.1 dex). The red squares and blue triangles represent stars with Jupiter-mass and Neptunian/super-Earth mass 
planets, respectively. The black circles refer to the stars without a planetary companion. The black dashed lines target the solar value. 



the NH stars have [Al/Fe], [Sc/Fe] and particularly [Mg/Fe] val- 
ues higher than 0, although JHs are also observed at much lower 
[X/Fe] values. We again note that the number of stars that exclu- 
sively host Neptunes and super-Earth planets is small, therefore, 
the conclusions regarding them should be considered with cau- 
tion. 



3.2. [X/Fe] versus [Fe/H] 

Previous planet-host studies (e.g. Neves et al. 2009; Petigura & 
Marcy 2011; Kang et al. 2011) have found that, in addition to 
iron, the abundances of various other metals are enhanced in 
these stars compared to stars with no known planets. Given the 
metal-rich nature of planet-hosting stars (higher overall iron con- 
tent), they are indeed expected to have a higher content of other 
metals as well. Some authors also referred to the existence of 
different [X/Fe] trends in planet-hosting stars compared to stars 
without planets for the same metallicity. For instance, potential 
differences were detected by Sadakane et al. (2002) for vana- 
dium and cobalt, by Bodaghee et al. (2003) for Ti, Mn, V, and 



Co, and by Gilli et al. (2006) for V, Co, Mg, and Al. Robinson 
et al. (2006) claimed to have found significant overabundances 
of [Si/Fe] and [Ni/Fe]. Gonzalez & Laws (2007) also found that 
Al and Si abundances are systematically lower for the planet- 
hosting stars in the higher metallicity region, and that the Ti 
abundance exhibits the opposite trend, which implies that the 
abundances of Na, Mg, Sc, and Ni might have some differences 
between planet hosts and stars without planets. Recently, Kang 
et al. (201 1), studying 34 PHSs and 18 comparison G-type stars, 
found that the [Mn/Fe] ratios of planet-hosting stars are higher 
than those of comparison stars over the entire metallicity range. 
These authors also found that [X/Fe] ratios of Mg, Al, Sc, Ti, 
V, and Co for PHSs are higher than those of comparison stars 
in metal-poor stars of [Fe/H] < -0.4 by more than 0.2 dex (they 
had only two stars with and two stars without planets in the men- 
tioned region). Most of these studies were in general limited to 
small samples of planet-hosting and stars without planets and as 
a result the general picture is so very difficult and contradictory. 

Knowing that different elements show different trends of 
abundances in the metal-rich and metal-poor domain in the 
Galactic chemical evolution, studing the [X/Fe] histograms is 
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Fig. 5. [ff/Fe] versus [Fe/H] (left) and [a/Fe] versus [Ref/Fe] (right) for the total sample. The red squares refer to the Jovian hosts and the blue 
triangles refer to the stars hosting exclusively Neptunians and super-Earths. The black dots represent the stars without planets. The black solid and 
green dashed lines represent the mean distributions of the planet host and non-host samples, respectively. The magenta dotted lines are the cubic 
fits for the low-cr (thin disc) and the high-a (thick disk + high-cr metal-rich) stars. The two metal-poor Jovian hosts in the red circles are the stars 
with low log g values discussed in Sect. 3.2. 



not sufficient to find small/subtle abundance differences between 
stars harboring planets with different masses and stars with- 
out planets. To explore the importance of individual elements 
in planet formation it is necessary to investigate the elemental 
abundances along with the metallicity. In this subsection, we will 
analyze whether there are any differences in the abundances of 
stars with and without planets for the same value of [Fe/H]. 

Fig. 4 shows the averages and standard deviations of [X/Fe] 
ratios for each [Fe/H] bin. The sizes of the bins are 0. 1 dex. The 
red squares and blue triangles represents stars with Jovian-mass 
and Neptunian-mass planets, respectively. The black circles re- 
fer to the stars without planetary companion. As shown in the 
figure, planet-hosting stars, both hosting Neptunes and Jupiters, 
are more abundant in a-elements than comparison stars when 
the metallicity is low. Giant-planet hosts have higher [Sc/Fe] 
and [Ti/Fe] values than stars without planets when [Fe/H] < - 
0.2 dex, higher [Al/Fe] and [Si/Fe] when [Fe/H] < -0.1 dex, 
and higher [Mg/Fe] ratios when [Fe/H] < 0. 1 dex. We note that 
the [X/Fe] values in the -0.5 < [Fe/H] < -0.4 dex metallicity 
bin for the Jupiter-hosting stars are not fitted with the general 
[X/Fe] vs [Fe/H] trends. In that bin we have three stars, two of 
which have "thin" disk [X/Fe] abundance ratios. Interestingly, 
these two stars are the most "evolved" planet hosts of our sample 
in that iron-poor regime. The stars, HD171028 and HD 190984, 
have 3.84 and 4.02 dex log g values, respectively, although the 
average log g of the metal-poor ([Fe/H] < -0.2 dex) planet hosts 
is 4.35 dex, and for the total planet-hosting sample it is 4.37 dex. 
It is difficult to conclude why these two stars show "peculiar" a- 
element abundances. Increasing the statistics will help to resolve 
this uncertainty. 

For low-mass planet host the increase in [X/Fe] ratios starts 
at higher metallicities than for high-mass planet hosts. Stars 
hosting Neptune-like planets have higher [Ti/Fe] values than 
their non-host counterparts when [Fe/H] < -0.2 dex, higher 



Table 2. The K-S probabilities that stars with planets of different mass 
and non-planet host stars have the same [X/Fe] distribution. Only stars 
with -0.2 < [Fe/H] < 0.4 dex are considered. 



Species 


Jovian hosts 


Neptunian hosts 




K-S Prob 


K-S Prob 


Na 


0.004 


0.242 


Mg 


0.010 


0.000 


Al 


0.009 


0.045 


Si 


0.012 


0.014 


Ca 


0.002 


0.328 


Sc 


0.049 


0.018 


Ti 


0.629 


0.763 


V 


0.123 


0.519 


Cr 


0.049 


0.099 


Mn 


0.002 


0.372 


Co 


0.031 


0.142 


Ni 


0.000 


0.195 


Fe 


0.000 


0.990 



[Al/Fe], [Si/Fe], [Co/Fe], [V/Fe] and [Sc/Fe] when [Fe/H] < 0.0 
dex, and they have higher [Mg/Fe] and [Na/Fe] (for Na the dif- 
ference is smaller) abundance ratios in all bins (except the most 
metal-poor and metal-rich bins where we have only one star with 
low [Na/Fe] ratio). Low-mass planet hosts also exhibit slightly 
higher [Ni/Fe] abundance ratios, although they are still consis- 
tent within the error bars in all metallicity bins except the -0.2 
< [Fe/H] < -0. 1 dex. With our relatively large and homogeneous 
sample we did not observe any abundance differences in man- 
ganese between stars with and without planets, as claimed in 
Kang et al. (2011), and we observed no nickel enrichment in 
stars with giant planets, as reported in Robinson et al. (2006). 

Based on these results, we can assume that some metals other 
than iron are involved in the process of planet formation, espe- 
cially when the amount of iron is lower than solar. Iron is not 
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the only abundant refractory element in the solar system. There 
are other fairly abundant elements (e.g. Mg, Si) with conden- 
sation temperatures comparable to iron (Lodders 2003; Lodders 
et al. 2009) that are very important contributors to the composi- 
tion of dust in planet-forming regions and represent the principal 
components of rocky-type planets. Thus, these results are con- 
sistent with the expectations, since a high [X/Fe] ratio means 
higher "global metallicity" and are moreover supported by the 
theoretical studies using the core-accretion model (e.g. Ida & 
Lin 2005; Mordasini et al. 2012). However, our analysis indi- 
cates that planet-hosting stars have higher [X/Fe] values than 
stars without planets for elements that are less abundant (e.g. 
Sc, Ti). This may indicate that not all solids are equally effective 
for planet formation. Alternatively, some elements (like Si or 
Mg) may "really" stimulate planet formation processes, although 
observed overabundance of otherelements may just reflect the 
similarities of the Galactic chemical evolution of the elements. 
Particularly, Si, Mg, Al, Sc, and Ti show similar [X/Fe] trends 
with [Fe/H], although in this case one might expect overabun- 
dance of Ca in planet host stars, which is not observed. More 
analyse are needed to understand which elements show a sec- 
ondary correlation with exoplanet presence. 

Although the observed abundance differences between the 
stars with and without planets are in general tiny compared to 
the standard deviations of the samples, they are systematic, and 
we believe that increasing the sample of planet-hosting stars will 
yield statistically more significant results. 



3.3. Exoplanets and the thick disc 

Haywood (2008, 2009) reported that at metallicities [Fe/H] < 
-0.3 most stars known to harbor giant planets belong to the 
thick disk rather than to the thin disc. Gonzalez (2009), using 
three samples of nearby stars hosting giant planets, showed that 
if compare stars with planets according to mass abundance of 
the refractory elements important for planet formation (instead 
of only iron), then thick disk and [Fe/H] -poor thin disk stars 
with planets have similar distributions. Studing this connection 
of planet-hosting stars with the thick disk is the subject of this 
subsection. 

In Fig. 5 we plot [a/Fe] against [Fe/H] and [a/Fe] against 
[Ref/Fe] for the total sample. The "Ref" index is the one pro- 
posed by Gonzalez (2009) and quantifies the mass abundances of 
Mg, Si and Fe. The red squares refer to the Jovian hosts and the 
blue triangles refer to the stars exclusively host Neptunians and 
super-Earths. The black dots represent the stars without planets. 
To illustrate the path of planet formation along the abundance 
ratios [a/Fe] we present the mean of the distribution (black solid 
line). For the comparison the mean of the [a/Fe] distribution of 
the stats without planets are also presented (green dashed line). 
The magenta dotted lines are the cubic polynomial fits for the 
low-c (thin disc) and the high-a stars (thick disk + high-a metal- 
rich, see Adibekyan et al. (201 1) for the chemical separation of 
the sample). 

With this data we confirm the finding by Haywood (2008, 
2009) that giant planet incidence is greater among the thick-disk 
population than among the thin-disk population for [Fe/H] < -0.3 
dex. We also show that this can be extended to the Neptunian 
hosts. At the mentioned metallicity region we have only three 
stars that harbor a planet in the thin disk (two of them are the 
aforementioned planet hosts with low log g values) and eight 
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Fig. 6. The same as Fig. 5 (left panel), but for stars with T tff =T Q ±300 
K. 



planet hosts in the thick disc 3 . Interestingly, this thin/thick pro- 
portion changes dramatically when we apply a purely kinematic 
approach to separate the thin and thick discs. The kinematical 
separation suggests that six or three host are from the thick disc, 
four or six stars belong to thin disk and one or two stars can be 
classified as transition stars depending on the kinematic criteria 
used - Bensby et al (2003) or Reddy et al. (2006), respectively 
(see also Adibekyan et al. 201 1). Moreover, following Gonzalez 
2009), and considering the Ref index instead of [Fe/H], we can 
see that the observed high abundance of metal-poor PHSs in the 
thick disk "disappears": most metal-poor planet-hosting stars in 
the thick disk have the same [Ref/H] distribution as their thin- 
disk counterparts (see the right panel of Fig. 5). This indicates 
that the thick-disk planet hosts might have less Fe than thin-disk 
hosts but they have relatively more Mg and Si, so they are not as 
metal -poor in terms of their ability to form planets. 

In the Fe-poor regime it is difficult to conclude abouts the 
main reason for most of the planet-hosting stars to lie in the high- 
a/thick-disk region. Generally, the Galactic thick disk is com- 
posed of relatively old stars (e.g. Bensby et al. 2005; Adibekyan 
et al. 201 1) that move in Galactic orbits with a large-scale height 
and long-scale length (e.g. Robin et al. 1996; Juric et al. 2008). 
At the same time, the thick-disc stars are known to have higher 
[a/Fe] ratios than the stars from thin disk at the same [Fe/H]. 
Below we present two indications that impel us to suppose that 
the certain chemical composition is the decisive factor, not the 
Galactic birth radius. 

• Fig. 5 clearly shows that planet-hosting stars show a con- 
tinuous increase in [a/Fe] with decreasing [Fe/H] at metallicities 
from -0.2 to -0.3 dex (starting from the thin disc, they rise to the 
thick disc) while the thin and thick disk stars without planets are 
separated very well by their [a/Fe] ratios (see Adibekyan et al. 
201 1). Fig. 6 better illustrates this separation; we plot only stars 
with effective temperatures close to the Sun by +300 K. In the 
metalicity region from -0.4 to -0.2 dex there is a clear separation 
between the two Galactic discs and there is no "transition" popu- 
lation between them, while the rise of [a/Fe] for planet hosts still 
exists. This fact may be an indication that planet-hosting stars 
follow a certain chemical evolution trend that requires a certain 
chemical composition for their formation. This does not mean 



3 The separations between the Galactic thin and thick discs were 
taken from Adibekyan et al. (201 1). 
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Fig. 7. Toomre diagram for the planethosting and stars without plan- 
ets. The red squares and blue triangles represent stars with Jupiter- and 
Neptune-mass planets, respectively. The magenta circles and green as- 
terisks refer to the high-a metal-poor (chemically defined thick disc) 
and high-o- metal-rich stars without planets, and the black dots refer to 
the chemically defined thin disk non-host stars. Dotted lines indicate 
constant peculiar space velocities, v pec = (U 2 LSR + V^ SR + W 2 SR ) i/2 = 20 
and 100 km s~\ 



that the planet-hosting stars constitute a separate stellar popu- 
lation, and only suggests that perhaps a planet can be formed 
everywhere the chemical requirements are satisfied. 

• The second hint, which may also indicate the importance 
of the a elements, is that planet-hosting stars start to have high 
[a/Fe] ratios at lower metallicites when they still belong to the 
thin disc. This is illustrated in the left panel of Fig. 5 where 
one can see that the average [oyFe] ratio for planet-hosting stars 
(black line) is higher than that for stars without planets (green 
dotted line) when [Fe/H] < - 0. 1 dex. 

In this subsection we considered low-mass and high-mass 
planet hosts together because, in general, they show similar 
[a/Fe] trends with the metallicity (although at "higher" metal- 
licities most of the Neptunian hosts lie above the "average" dis- 
tribution of the combined host sample). 

The fact that most of the metal-poor planet hosts belong 
to the thick disk leads to an interesting conclusion. This sug- 
gests that when we calculate the frequency of planet hosts at low 
metallicities we should consider only stars belonging to the same 
Galactic population. In particular, the frequency of planet hosts 
in the metallicity region from -0.6 to -0.3 dex is about 5%, but 
when we consider only stars belonging to the thick disc, we will 
have about 14% and in the thin disk about 2%. To compare, the 
frequency of planet hosts in the metallicity region from -0.3 to 
dex in the thin disk is about 9%. This may also suggest that the 
observed paucity of very metal -poor planet hosts just reflects the 
low density of the thick disk stars in the solar neighborhood. 



Table 3. The average values of the Galactic orbital parameters and ages 
for the stars with giant planets, stars that exclusively host Neptunians, 
stars without planets, and stars without planets with \ovi-v pec . The stan- 
dard errors of the mean values are also presented. 
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Ecc 


Ron/ 


Age 


NHs 


0.33±0.06 


0.19±0.02 


8.00±0.001 


5.5±0.3 


JHs 


0.21±0.02 


0.14±0.01 


7.99±0.002 


4.8±0.3 


Non-hosts 


0.29±0.02 


0.15±0.00 


7.99±0.001 


5.2±0.1 


\cw-v pec 


0.10±0.01 


0.05±0.00 


7.98±0.0()3 


4.8±0.3 



4. Kinematics of planet-hosting stars 

So far, more than a dozen studies were published focused on the 
kinematics of the PHSs and their relation to different stellar pop- 
ulations. The first papers did not find any significant kinematic 
peculiarity of PHSs (e.g. Gonzalez 1999; Reid 2002; Barbieri 
& Gratton 2002). Conversely, Haywood (2008;2009), combin- 
ing the chemical and kinematic properties of the PHSs, con- 
cluded that most metal-rich stars that host giant planets originate 
from the inner Galactic disk. The same scenario for the origin of 
metal-rich planet hosts were proposed in some other papers (e.g. 
Ecuvillon et al. 2007; Santos et al. 2008; Viana Almeida et al. 
2009). 

To study the kinematics of planet-hosting stars in Fig. 7 
we presented their positions in the UW vs V velocities space 
(Toomre diagram). For details of the computation of the space 
velocity components we refer the reader to Adibekyan et al. 
(2012, submitted). In the diagram we used only stars without 
planets with [Fe/H] > -0.6 dex, because we found no planet be- 
low this value. Dotted lines indicate constant peculiar space ve- 
locities, v pec = (U 2 LSR +Vl SR +Wl SR ) l/2 = 20 and 100km s _1 . We 
first note is that there is no kinematically very "cool" Neptunian 
host with Vp ec < 20 km s~ l , while 18 stars that host Jovian plan- 
ets (~17 %) have peculiar space velocities slower than 20 km 
The number of stars without planets and v pec < 20 km s~ l 
is 121, which is »13 % of the total sample of non-hosts ([Fe/H] 
> -0.6 dex). 

First we checked if the different groups lie at the same dis- 
tances from the Sun. We found that NHs are, on average, closer 
to the Sun than JHs or stars without any planet companions. The 
average distance of NHs is »18 pc, while the average distance of 
JHs is «32 pc, for stars without planets it is «38 pc, and for the 
stars with v pec < 20 km s~ l is ^39 pc. This result is expected, 
because low-mass planets are easier to find at smaller distances 
due to the higher apparent magnitudes of their hosts. However, 
we have found that the difference in the distances is not the rea- 
son for the "high" peculiar velocities of NHs. Among the stars 
without planets that are closest to the Sun (the average distance 
of this subsample was x;18 pc) we found 18 stars (^16% of the 



the total subsample) with v pec < 20 km s . 

To understand the root of this kinematical "peculiarity" of 
NH stars we collected the main orbital and kinematical proper- 
ties of these stars and compared them to those from other groups 
of stars. To obtain the Galactic orbital parameters of the stars, we 
cross-matched our sample with the Geneva-Copenhagen Survey 
(GCS) sample (Casagrande et al. 2011), which provides the ec- 
centricities of the orbits, the maximum vertical distance (Z max ) a 
star can reach above the Galactic plane, the Galactic radial po- 
sitions, and the ages of about 700 of our stars (23 NHs, 77 JHs, 
548 non-host stars with [Fe/H] > -0.6 dex, and 55 non-hosts with 
[Fe/H] < -0.6 dex). The distributions of the Galactic orbital pa- 
rameters and space velocities of the stars with low-mass plan- 
ets (blue shade), with high-mass planets (red line), stars without 
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Fig. 8. Distributions of the Galactic orbital parameters and space velocities of the stars with Neptune-mass planets (shaded blue), with Jupiter-mass 
planets (red dashed line), stars without planets and [Fe/H] > -0.6 dex (black dotted line), and stars without planets with [Fe/H] > -0.6 dex and 
v pec < 20 km s~ l (green dotted-dashed line). The distributions of different groups were set lower/higher for the sake of clarity. 



planets with [Fe/H] > -0.6 dex (black line), and the non-host 
stars with [Fe/H] > -0.6 dex and v pec < 20 km s~' (green dot- 
ted line, hereafter low-v pec stars) are presented in Fig. 8. The 
average values of the physical parameters aand the ages for the 
groups of stars are presented in Table 3. 

Inspection of the figure and the associated table shows that 
NH stars, on average, have higher Z max , higher eccentricities, 
and higher Re,,/ than their high-mass planet-hosting counter- 
parts and stars without planets. The differences, as expected, 
are higher when compared to the low-y ; , t , c stars. From the bot- 
tom panels of Fig. 8 we can see that NHs have lower Vlsr and 
Wlsr space velocity components than stars from other groups. 
All these differences indicate that NHs tend to belong to the 
"thicker" disc, although JHs show the same frequency in the 
"thinner" disk as the stars without planets in our sample. This 
finding confirms also the last column of Table 3 where one can 
see that NH stars are on average older than the JHs (although the 
difference is small). 

It is difficult to conclude about the main reason of the ob- 
served tendency of NH stars to belong to the "thicker" disc. In 
Fig. 7 we can see that NHs follow the distribution of high-cr stars 
in the UW vs V velocity space. To see if the NHs and stars with- 
out planets with the same metallicity (-0.6 < [Fe/H] < 0.25) and 
kinematics (v pec > 20 km s~ l ) have the same abundances of a- 
elements, we compared their [Mg/Fe], [Si/Fe] and [Ti/Fe] ratios. 
We found that the NH stars are more enhanced by Mg (the differ- 
ence in [Mg/Fe] * 0.06 dex) than by Si (A[Si/Fe] « 0.02 dex) and 
Ti (A[Ti/Fe] « 0.02 dex), and both groups of stars are enhanced 
by these elements compared to the low-v pec stars without plan- 
ets. This may be a subtle hint confirming our previous finding 
(see Sect. 3) that to form Neptune-mass planets, some elements, 
like Mg, are particularly important. 



5. Conclusions 

We presented a differential abundance analysis between planet- 
hosting and stars without planets for 12 refractory elements (Na, 
Mg, Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) for a total of 
1111 nearby FGK dwarf stars. Of these stars, 109 are known to 
harbor high-mass planetary companions and 26 stars are exclu- 
sively hosting Neptunians and super-Earths. The precise spectro- 
scopic parameters for the entire sample were taken from Sous a 
et al. (2008, 201 la, 201 lb) and the abundances were taken from 
Adibekyan et al. (2012, submitted). 

The inspection of the [X/Fe] histograms suggests that in gen- 
eral, the [X/Fe] distributions of NHs start at higher [X/Fe] values 
compared to the distributions of giant-planet hosts. This shift to- 
ward the higher [X/Fe]values is most clearly observed for Mg. 
For Co, Na, Ni, V, and Mn we observe that there is a general 
increase in the frequency of stars with giant planets, with in- 
creasing [X/Fe]. We also found that most of the Neptunian hosts 
have [Al/Fe], [Sc/Fe] and [Mg/Fe] values higher than 0, although 
Jovian hosts can have much lower [X/Fe] values for these ele- 
ments. 

With our relatively large and homogeneous sample of planet- 
hosting and stars without planets we found that the [X/Fe] ra- 
tios of Mg, Al, Si, Sc, and Ti for giant hosts are systematically 
higher than those of comparison stars at [Fe/H] < -0.1 +0.1 
dex. Simultaneously, stars hosting Neptune-like planets show 
higher [Ti/Fe] ([Fe/H] < -0.2 dex), [Si/Fe], [Al/Fe], [Co/Fe], 
[V/Fe] and [Sc/Fe] ([Fe/H] < 0.0 dex), and higher [Mg/Fe] and 
[Na/Fe] (over the entire metallicity range) values than their non- 
host counterparts. Low-mass planet hosts also exhibit slightly 
higher [Ni/Fe] abundance ratios, although they are still consis- 
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tent within the error bars in all metallicity bins except the -0.2 < 
[Fe/H] < -0.1 dex. 

We confirmed that planet incidence is greater among the 
thick-disk population than among the thin disk for [Fe/H] < - 
0.3 dex. At lowest metallicities we observed only two stars with 
low [ayFe] ratios and both have lower log g values (3.84 and 
4.02 dex). 

We also discussed recent debates about the high abun- 
dance of metal-poor planet hosts in the thick disk (Haywood 
2008, 2009; Gonzalez 2009), i.e. the main reason that most Fe- 
poor stars lie in the high-cr/thick-disk region - is it a special birth 
place in the Galaxy or a certain chemical composition. Our re- 
sults allow us to propose that the certain chemical composition 
and not the Galactic birth place is the determining factor. 

The study of kinematical properties of planet host stars 
shows that there is no kinematically very "cool" low-mass planet 
host star with a peculiar velocity slower than 20 km s -1 , while 
~17 % of stars that host high-mass planets have v pec < 20 km 
s . Although we found that NHs are closer to the Sun than JHs 
and stars without planets, this is not a reason of the lower veloc- 
ity limit. Inspecting of the Galactic orbital parameters of these 
stars shows that NHs tend to belong to the "thicker" disk com- 
pared to their high-mass planet-hosting counterparts. We also 
found that they follow the distribution of high-cr stars in the UW 
vs V velocity space, but they are more enhanced in Mg than 
high-cr stars without planets. Summarizing the results obtained 
for the Neptunian-like planet hosts, we can assume that to form 
low-mass planets some elements, like Mg, might be particularly 
important. 
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